Weizsacker energy of unitary Fermi gas in a harmonic trap 

Alexander L. Zubare^^*] 
Department of Physics, Purdue University 
West Lafayette, Indiana 47907 

The universal method of construction of the rigorous lower bounds to the Weizsacker 
energy is presented. We study a few-fermion system at the unitarity Upper and 
lower bounds to the density functional theory (DFT) ground state energy within 
the local density approximation (LDA) are given. The rigorous lower bounds to 
the accuracy of the method are derived. 



PACS numbers: 31.15.Ew, 71.15.Mb, 03.75.Ss 



* e-mail: zubareva@physics.purdue.edu 



There has been a lot of interest in systems of fermions at the unitarity [1-3] (when the 
scattering length diverges, the Bertsch many-body problem [4]). While Refs.[5-8] consider 
homogeneous systems, Refs.[9-13] present ab initio calculations of the properties of trapped 
fermionic atoms. 

The modern DFT is based on the Kohn-Sham approach [14], where the noninteracting 
kinetic energy, T, is calculated in terms of the Kohn-Sham orbitals, although Ref. [15] proved 
the basic existence of functional T(p), where p is the density, / p{f)d^r = N, and is 
the particle number. The accurate density-functional approximation to the kinetic energy 
would reduce dramatically complexity of the DFT calculations (here we note the superfluid 
extension of the DFT given in Refs. [16-18]). For applications of the DFT to the nuclear 
structure physics see web site, constructed for the universal nuclear energy density functional 



(UNEDF) collaboration, http://unedf.org 



The kinetic energy functional cam be written as 



T[p] = ^ I TW))d\ (1) 



2m 

where the semiclassical expansion for the kinetic energy is [19-23] 

^2 



T[p\ = —j {tMp) + r,{p) + n{p) + ...)d\ (2) 



where 



tMp) = \{2,^-?l'p'l\ (3) 
r.M = i<?^, (4) 

^'^P^ = ^(37r2)"^V/^[8(^)^ - 27(^)^^ + 24(^)2]. (5) 
o48U p P P P 

The Weizsacker energy T]y[p] [24] 

is an important component of the DFT kinetic energy. Indeed, Tw is considered exact in 
the limit of rapidly varying density p [25,26] and nine times large than the second term (the 
Kirgnitz correction [19]) of the semiclassical expansion (2), which have to be considered as 
an asymptotic expansion [23]. 



We present in this letter the universal method to construct the rigorous lower bounds to 
the Tw which allows us to study a few-fermion system at unitarity. 
We rewrite Eq.(6) as 

Tw[p] = / [|^[V/^'(0]' + (V(f, Ai, A2, ...Am) - V(f, X,, A2, ...AM))]p(OdV, (7) 

where V{f, Ai, A2, ...Am) is an one-particle potential and A are parameters. 
Introducing an auxiliary Hamiltonian Haux as 

Haux^-:z-A + V{f,X,,X2,...XM), (8) 
2m 

and using 

< ij\Haux\ip >> Eg < >, (9) 

we get 

Tw[p] > NEciXi, A2, ...Am) - / Vif, X,, X^, ...XM)p{^d\ (10) 

where Eg is the ground state energy of the auxiliary Hamiltonian (8). Therefore, a set of 
optimal values of parameters Aj which yield an optimal value for a lower bound to the 2V[p] 
is given by 

T'^^'^'ip] = max[NEG{X,, A2, ...Am) - / V{f, X,, A2, ...AM)p(r)dV] (11) 

Now let us collect some lower bounds to the Weizsacker energy Tiy[p]: 

(i) V{f, A) = mX^ry2, Eg{X) = {3/2)hX, A > 0, 

Tw[p] > 12 

(ii) V{f,X) = -A/r, EGiX) = X^m/{2h^), A > 0, 

(iii) V{f, Ai, A2, A3) = {m/2){Xix' + Xly' + X^z^), EGiX,, A2, A3) = ih/2){X, + A2 + A3), 
Xi > 0, 



Twip] > '^i^::^ + — ^ + — ^). (14) 



h'^N^ 1 1 1 

Sm < x"^ > <y^ > < > 
(iv) V{f, A) = mAVV2 + f3{(3 + l)ti^ / {2mr^) , Eg = nX{3/2 + f3), (3 > 0, X > 0, 



^ r 1 ^ 2 3iV2+ < >< > 

TV P > 7^ < ^ > 77 o 15) 



(v) V(r, A) = -A/r + P{P + l)Tf l(2mr''), Eo = -\^m/(21i'(l + Pf), /J > 0, A > 0, 

N<r-^>> 

where < >= / p{r)r°'d^r. 

The bounds (12), (13), (15) and (16) where previously proved in Refs. [27,28], while the 
bound (14) is new. 

In Refs. [29,30] and later in Refs. [31-42] the dynamics of strongly interacting dilute Fermi 
gases (dilute in the sense that the range of interatomic potential is small compared with 
inter-particle spacing) consisting of a 50-50 mixture of two different states and confined in 
a harmonic trap Vext{^) = {m/2){uj\{x'^ + y'^) + ujIz^) is investigated in the single equation 
approach to the time-dependent density-functional theory. 

Let us come back to the variational formulation of the Kohn-Sham time-dependent theory 

5 j dt< ijjlindt - H\ijj >= 0, 17) 

where 1-0 > is a product of two Slater determinants, one for each internal state built up by 
the Kohn-Sham orbitals ipi, ^-nd H — T + U is the LDA Hamiltonian. 
Using two approximations 

(i) local transform ipi Ki 0j exp(i^x/m), where ipi and x real functions, 
and 

{n)<(l>m>^SlI^TF{p)d'r + Tw[p], 

where |0 > is the product of two Slater determinants built on 0j alone, we can derive the 
DFT equation of Ref.[29] 

th— = -— V^* + Kxt* + (18) 
at 2m 

where Vxc{f,t) — [ ^^gp^^ ]p=p(f,f)) ^(p) is the ground-state energy per particle of the homo- 
geneous system and p{f,t) —\ ^{f,t) p . In the case of large but finite number of atoms 
N, at small distances the ratio | Vp | / p^^^ is small and both the Kirzhnitz correction and 
the Weizsacker correction are negligible. On the contrary, the Weizsacker correction is ex- 
pected to determine the asymptotic behavior of the density at large distances. As for the 
case of relatively small number of atoms, we expect that the Kirzhnitz correction would be 
a reasonable approximation to the kinetic energy 

T[p\^TTF[p\ + \Tw[p\d\ (19) 



= + 2^-**^ + 214c*., (20) 



where rrf[p] = £/(TTF(p)rfV. 

Recently, Ref.[41] has considered the following nonhnear equation 

t 

dt Am 

with = p/2. 

For the stationary case Eq.(20) reduces to 

_ V^* + K.t* + Kc* = (21) 
8m 

which corresponds to the following approximation of the kinetic energy 

T[p\ ^ Ttf[p\ + -Tw[p] (22) 



with K — 4. 

For the remainder of this letter we will test the accuracy of the approximation (22) with 
K = 9 and k = 4 for few-fermion systems at unitarity in a spherical harmonic trap 



2 2 

K.t(r) = . (23) 



In the approximation (22), the ground state energy is given by the minimum of the energy 
functional 

.]{% = ^Tw[p] + I Ve^tpd'r + I e{p)pd\ (24) 

where e(p) = (1 + P)3Ji? kp / (lOm) , kp — (Stt^p)^/^, p — and the universal parameter /3 
is estimated to be /3 = —0.56 [7]. 

Introducing an auxiliary Hamiltonian 



we can rewrite Eq.(24) as 



= ^Tw\p\ + / (Kxt - ^)pc?V + j{H + e{p))pd^r. 



Omission of {1/k,)Tw[p] + J{Vext — H)pd?r yields our approximation for the ground state 
energy 

E ^Nhu;\ - + Etf{VI - Xuj), (26) 
2 V 
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where Etf is the Thomas- Fermi energy which is given by Etf{^) — + Phu;{3N)^/^/4: 

Projecting |^ >on the complete basis states \n >, obtained from h\n >— en\n >, where 
h = -{h'^/{2Km))A + we get 

3 /a 

< >= < >< >> -NTujj\1-. 

n 2 \ K 

Therefore, we conclude that our approximation for energy, given by Eq.(26), is a lower bound 
to the ground sate energy, Eq.(24). The optimal value of parameter A which maximizes the 
energy, Eq.(26), will yield an optimal value of the lower bounds to the ground state energy 
given by 

h^-^-27^V^ 4 ■ ^^^^ 

For large N, the finite N correction can be written as 

E(-) _ {3Nr/^VlTl3 (3iV)V3 

hoj 4 2k^/TT^ ^ ^ 

Bhaduri, Murthy and Brack [43] have recently presented a semiclassical approximation, 
assuming the particles obey the Haldane-Wu fractional exclusion statistics at unitarity 

where the statistical parameter g is related to the universal parameter /3 by the relation 
g — {1 + We note here that for the case of k = 4 [41] our lower bound, Eq.(28), agrees 

with the expression of Ref . [43] . 

To calculate upper bounds, E^~^\ we employing Fetter's trial functions [44] 

p'/\r}^c[l-(l-q)M']Tk, (30) 

where 7, q are the variational parameters and c is the normalization constant, to minimize 
the functional J, Eq.(24). 

From Table I, we can see that E'^-^ > E^*^ for k = 4 and 2 < N/2 < 14. Therefore 

\E{k = A)-E^'^\ 

EMC ^>A(/. = 4), 

where A(fi: = 4) = (£"("'(«; = 4) - E^^)/E^'^ and E{k) is the exact ground state solution 
of Eq.(24). For K = 9 and 3 < N/2 < 14, < E^^ (see Table I), that is why 

\E^^ -E(k^9)\ 

E^ ^>A(«: = 9), 

6 



where A{k = 9) = {E^^ - E^+\k = 9))/E^^. 

We can, therefore, state that A(«;) is the rigorous lower bound to the accuracy of the 
approximation (22) with k = 4, 2 < N/2 < 14 and k = 9, 3 < N/2 < 14. As for the lower 
and the upper bounds, they provide the actual solution of Eq.(24), + E^^y2 

within ±6 accuracy, with 6 < 1% for both k, = 9 and n = A and for 2 < N/2 < 15. 

The predictions of Eq.(24) for the ground state energy and results of Ref. [10] are shown 
in Fig.l. A very good agreement between the Kirgnitz approximation, k — 9, and the ab 
initio calculations of Ref.[10] can be seen for 2 < A^/2 < 4. However, for 12 < A^/2 < 15 the 
K — A approximation gives better results than the k — 9 approximation. 

Finally, we note that the kinetic energy functional assumed here are not unique. In our 
future work we will consider other possible forms. 

In summary, we have constructed the rigorous lower bounds to the Weizsacker energy. 
As example of application, we have studied few-fermion systems at unitarity consisting of 
50-50 mixture of two different states and confined in a spherical harmonic trap. The rigorous 
lower bounds to the accuracy of the method are derived. We have tested the kinetic energy 
functionals by comparisons with ab initio calculations and have found that while the second 
order gradient expansion is a very accurate for relatively small N, Fig.l indicates that the 
K — A approximation provides significantly better results for larger N. 
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Table I. The energies \ E'^'^^ and the energy calculated within the fixed- node diffusion 



Monte Carlo method, E^^ [10], all 


in units of fiu; for A?" < 30 [see 


the text for further details] . 
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Fig.l. Ground state energy per particle of few-fermion systems at unitarity in a spherical 
harmonic trap in units of Tiuj as a function of number of atoms N . The dashed and the 
dotted-dashed lines represent {E^~^ + (£'(+))/2 for k = 9 and k = 4, respectively. The 
circular dots indicate results of Ref.[10]. 
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